Quasiparticle bands in cuprates by quantum chemical methods: 
towards an ab initio description of strong electron correlations 
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Realistic electronic-structure calculations for correlated Mott insulators are notoriously hard. 
Here we present an ab initio mu/ficonfiguration scheme that adequately describes strong correlation 
effects involving Cu 3d and O 2p electrons in layered cuprates. In particular, the O 2p states giving 
rise to the Zhang-Rice band are explicitly considered. Renormalization effects due to nonlocal spin 
interactions are also treated consistently. We show that the dispersion of the lowest band observed 
in photoemission is reproduced with quantitative accuracy. Additionally, the evolution of the Fermi 
surface with doping follows directly from our ah initio data. Our results thus open a new avenue 
for the first-principles investigation of the electronic structure of correlated Mott insulators. 
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Angle-resolved photoemission spectroscopy (ARPES) 
has proved to be an invaluable tool for probing the cor- 
related electronic structure and the one-particle spec- 
tral function, A(k, w), of interacting many-body systems 
such as the Sd transition-metal (TM) oxides. In addition 
to the ubiquitous strong d-d Hubbard correlations, for 
late TM oxide compounds like the cuprates and nicke- 
lates, the electron-removal spectra are strongly affected 
by 2p- TM Sd charge transfer and rehybridization ef- 
fects. Detailed analysis of the photoemission data, see 
Ref. [H , led to the conclusion that the valence-band states 
at low binding energies actually have predominant O 2p 
character in nickelates and cuprates. Of particular inter- 
est is the nature of the first electron-removal states in the 
high-Tc superconductors. For the (Mott) insulating par- 
ent compounds, the formal valence states of copper and 
oxygen are Cu^+ Sd^ and 2p^, respectively, with 

one 3d hole at each Cu site. In a seminal work, Zhang 
and Rice (ZR) [3| argued that extra holes doped into the 
Cu02 layers mainly populate the ligand 2p levels and give 
rise to impurity-like states involving singlet coupling with 
Cu 3d holes. This was indeed confirmed by both quan- 
tum chemical wave-function based calculations on finite 
clusters 0, 01 and periodic density- functional calculations 
within the LDA-|-[/ approximation |5|]. Matrix elements 
associated with the hopping of the 2p hole to neighboring 
plaquettes were also computed in Ref. 0]. However, the 
moderate size of the clusters (up to four CUO4 plaquettes 

) precluded a study of the effect of the antiferromag- 
netic (AF) spin background on the effective hoppings. 

Calculations based on t-J model Hamiltonians pre- 
dict that the dispersion of the quasiparticle bands is 
strongly renormalized by the underlying AF interactions 
to 0(J) 0. A relatively small width for the bands 
close to the Fermi level was indeed observed by ARPES 
investigations. ARPES thus plays a central role in study- 
ing the dressing of charge carriers by interactions with the 
spin (and lattice) degrees of freedom 0, [1| • In addition, 
these measurements provide detailed information on the 
doping dependence of the Fermi surface (FS). In par- 
ticular, the "nodal- antinodal" dichotomy in underdoped 



cuprates 0] was also emphasized by very recent Raman 
scattering data the intensity along the nodal di- 

rection scales with Tc, while the antinodal signal actu- 
ally anti-correlates with Tc(x). The details of the FS are 
important for the proper interpretation of other exper- 
iments as well, such as the inelastic neutron scattering 
, and put strong constraints on theoretical models for 
cuprates. 

To date, the above issues have been systematically ad- 
dressed by cluster dynamical mean-field investigations 
on effective one-band Hubbard models [3]- 
A first-principles correlated electronic-structure calcula- 
tion, however, is still quite a distant goal, especially for 
late TM oxides. Here, we present an ah initio quantum 
chemical study, describing how the detailed quasiparticle 
dispersion can be obtained with quantitative accuracy. 
To this end, multiconfiguration (MC) calculations were 
performed on clusters that are large enough to account 
for spin interactions in the neighborhood of the ZR hole. 
Renormalization effects on the nearest-neighbor (NN) ef- 
fective hopping were previously investigated at the ab 
initio level in Ref. In this Letter, we show that ac- 
curate estimates for the longer-range hoppings are crucial 
for explaining both, the dispersion observed by ARPES 
measurements, as well as the Fermi surface, in detail. 

The ab initio calculations are performed on finite frag- 
ments cut from the periodic system. In constructing the 
MC wave-function, the orbitals of each finite cluster are 
partitioned into three different sets: the inactive levels, 
doubly occupied in all configurations, the virtual orbitals, 
empty in all configurations, and the active orbital set, 
where no occupation restrictions are imposed. This type 
of MC wave-function is also referred to as a complete- 
active-space (CAS) wave-function [Tg]. In the case of an 
undoped cluster, for example, with formal Cu 3d^ and 
O 2p^ occupations of the valence levels, the most natural 
choice for the active space is to include all copper ?>d,^2_y2 
orbitals, akin to the wave- function in the one-band Hub- 
bard model at half filling. Still, all integrals, including 
the on-site and longer range Coulomb interactions, are 
here computed ab initio. Although electrons in lower 
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(inactive) energy levels, e.g., O 2s and 2p, do not par- 
ticipate explicitly in the construction of the many-body 
expansion of the MC wave-function, they are still allowed 
to readjust to fluctuations within the active orbital space. 
If extra holes are created, the active space must be en- 
larged with orbitals from the inactive orbital group. A 
single doped hole requires one orbital to be transferred 
from the inactive to the active space. Since different elec- 
tronic states have in many cases very differently shaped 
orbitals, it is important to express each state in terms 
of its own, separately optimized set of orbitals. For the 
lowest electron-removal state, for example, the orbital 
added to the active space turns into a ZR p-d compos- 
ite in the variational calculation, localized on a given 
Cu04 plaquette. The character of the active orbitals may 
change substantially, however, when separate CAS self- 
consistent-field (CASSCF) optimizations are performed 
for higher-lying hole states, see below. 

The effective hoppings (and quasiparticle bands) are 
computed by using the overlap, Sij, and Hamiltonian, 
Hij, matrix elements between (A'^— 1) wave-functions hav- 
ing the extra hole located on different plaquettes 
of a given cluster. Each of these (A^— 1) wave-functions, 
is obtained by separate CASSCF optimizations. 
A similar approach was previously applied to simpler, 
noncontroversial systems such as MgO (see, e.g., 
and references therein). For degenerate {Hii = Hjj) hole 
states: t = 1/2(6^ - e,) = (Hy - S,jHu)/{l - Sf^), 
where and ej are the eigenvalues of the 2x2 secu- 
lar problem. For non-degenerate (Hu^Hjj) hole states, 
t = l/2[(ej - ei)^ - {Hjj - Huf]^/^. The overlap and 
Hamiltonian matrix elements between the individually 
optimized, nonorthogonal CASSCF wave-functions are 
obtained by State-Interaction (SI) [l^ calculations. 

The short-range magnetic correlations are described by 
including in our finite clusters extra CUO4 units around 
those plaquettes involved in the hopping process. We call 
those latter plaquettes "active" . All-electron basis sets of 




TABLE I: Occupation numbers (ON's) of the Cu d^2_y2, Cu 
d^^2_^2, in-plane cr-type O Px/py, and apical (ap.) O Pz 
atomic orbitals for the ZR and the d^z2_^2 hole states. For 
undoped plaquettes, the ON's of the a 2p orbitals are ~1.85. 
The "missing" charge is in the Cu 4s and 4p orbitals. 



ON's 




d^^2_r2 


cr Px/Py 


ap. Pz 


ZR hole 


1.05 


2.00 


1.60 


1.95 


d^^2_j.2 hole 


1.40 


1.15 


1.70 


1.85 



triple-zeta quality [l9| are applied for the Cu and O ions 
of the active plaquettes. The core electrons of the re- 
maining ions in each cluster are modeled by effective core 
potentials (ECP's) [23|. That this approximation works 
very well was shown for the NN hoppings in Ref . [l^ . To 
describe the finite charge distribution at the sites in the 
immediate neighborhood of the cluster, we model those 



ions by effective ionic potentials, see [21[. Beyond these 
neighbors, we use large arrays of point charges that repro- 
duce the Madelung field within the cluster region. The 
various clusters employed in our calculations are sketched 
in Fig. 1. Apical ligands are explicitly included in our 
calculations only for the active plaquettes. Other apex 
oxygens are represented by formal point charges. All cal- 
culations were performed with the molcas package [IJl . 
The structural data reported for the Lai.85Sro.i5Cu04 
compound [2;|| were used, with an in-plane lattice con- 
stant a = 3.78 A. 

The dispersion of d-like states on a square lattice is 
given by the following relation: e(k) = —2t (cos k^a + 
cos ky a)— ^t' cos kx a cos kyU — 2t" {cos 2kxa + cos 2ky a)... , 
where t, t' , t" are the hopping integrals between NN, 
second-NN, and third-NN sites and each effective site is 



FIG. 1: Sketch of the finite clusters employed for the calcula- 
tion of the effective hoppings of the ZR-like hole. 



a whole CUO4 plaquette. As discussed in Ref. [15[, a hole 
in an O 2p a-type orbital induces ferromagnetic (FM) 
correlations among the NN Cu 3d spins. These spin po- 
larization "tails" around each 2p hole might actually play 
a role in pairing. When moving through the AF back- 
ground, the hole must drag along this spin polarization 
cloud at nearby Cu sites, which gives rise to a substan- 
tial reduction of the effective hopping parameters. The 
spin polarization and relaxation effects at nearby sites 
are treated consistently in our approach because a sepa- 
rate optimization is performed for each particular |\E'f 
configuration. A renormalization by a factor of four was 
found for the NN effective hopping, from = 0-5— 0.6 eV 
to t = 0.135 eV [15]. The renormalized values for t' and 
t" , obtained by CASSCF and SI calculations on clusters 
like those shown in Fig. 1, are —0.015 and 0.073 eV, re- 
spectively. Surprisingly, t" turns out to be about half 
of the renormalized NN hopping t. The effect of such a 
large third-NN hopping is illustrated in the upper panel 
of Fig. 2. It leads to a weakly dispersing band at the 
(7r,0) point, with a small "dip" around that region, and a 
local energy maximum not far from the (7r/2,7r/2) point. 
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For comparison, we also plotted in Fig. 2 the dispersion 
along the nodal direction for t" = (the dotted curve) . 

An important issue that must be clarified is whether 
there is significant mixing between the lowest-energy hole 
state, the ZR-like state, and hole states at higher bind- 
ing energies. LDA+[/ calculations [H| indicate that the 
separation between the ZR configuration and the next 
electron-removal state is of the order of 0.1 eV. This sec- 
ond state is related to the creation of a 3^322 hole. Us- 
ing model Hamiltonian calculations, Eskes and Sawatzky 
[26| predicted a separation of about 1.5 eV. They also 
pointed out that the largest mixing should occur along 
the antinodal (0,0)-(7r,0) direction. 

Quantum chemical calculations Q were previo usly 
used for the interpretation of the optical absorption [27| . 
The shoulder at 1.5 — 2.0 eV was assigned to transitions 
between the d^2_y2 and d^^2_^2 levels |27|. Regarding the 
electron-removal spectrum, our CASSCF calculations in- 
dicate a separation Ae = 0.60 eV between the ZR and the 
d3z2_r2 (A^ — 1) states. Occupation numbers of the rele- 
vant atomic orbitals (MuUiken charges, see [l3|) are listed 
in Table L For the (13^2 hole state, the two holes on 
a given plaquette are high-spin coupled. Due to the mu- 
tual Coulomb repulsion and strong p-d hybridization, the 
{x^—y^) hole has large weight at the ligand sites. The NN 
d-d spin correlations are AF for this higher excited state. 
Highly accurate results for the on-site energies of the ZR 
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FIG. 3: Three-dimensional plot of the ZR-like quasiparti- 
cle dispersion, including the effect of the ZR-z^ interaction. 
Constant-energy contours are also drawn in the figure. The 
dashed line resembles the "FS" measured in doped La2Cu04 
[25! and other cuprates 0, see the sketch in the inset. 



larger than 20% in TM oxides. Illustrative in this re- 
spect are the results of Martin and Hay [3| for the d-d 
excitations in undoped clusters. For the inter-site ma- 
trix elements, the corrections to the CASSCF results are 
much smaller, see. 
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e.g., the discussion in Refs. 
We found that the mixing between ZR and hole 
states on the same plaquette is negligible. However, the 
inter-site matrix elements are large. The SI calculations 
indicate a NN effective hopping of 0.15 eV between the 
ZR and (A^— 1) states. We denote this quantity as tm- 



and d'iz'^_r2 hole states require calculations beyond the The k-dependent term describing the mixing between the 



CASSCF level. It is known, however, that the correc- 
tions brought by more sophisticated techniques such as 
the multi-reference configuration-interaction method [l6l | 
to the on-site CASSCF relative energies are usually not 
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FIG. 2: (upper panel:) The lowest electron-removal band, 
without including the interaction with the d2,^2_^2 hole state 
(dashed line) and after including this interaction (solid line). 
There is no mixing along the nodal direction. Close to the 
r point, the two bands cross. For simplicity, only the lowest 
electron-removal band is shown. The reference is the value of 
the on-site Hamiltonian matrix element of the ZR state, see 
text. The dotted line along the nodal direction corresponds 
to t" = 0. (lower panel:) ARPES data obtained for La2Cu04 
by Ino et al. [Ij. Units of eV are used in both plots. 



two bands, ZR and z , now reads tm{coskxa — coskyO). 
The NN hopping between degenerate hole states is 
much smaller, 1^2 = 0.01 eV. It is now trivial to diago- 
nalize the k-dependent 2x2 matrix to yield the renor- 
malized quasiparticle bands. The resulting dispersion of 
the "(a;^ — y^)" band is compared to the ARPES data in 
Fig. 2. It is known that as the F point is approached, 
spectral weight is transferred from the lowest binding- 
energy peak to a higher-energy, rapidly dispersing feature 
whose origin is a matter of active debate. This higher- 
energy feature, not visible in the lower panel of Fig. 2 
but clearly observed in Ref. [2^ and other studies, and 
the spectral-weight transfer at the F point are not ad- 
dressed in this work. Remarkably, in all other respects, 
our ab initio data agree very well with the dispersion 
of the lowest ARPES band. In particular, the flat dis- 
persion around (7r,0), the maximum near (77/2, 7r/2), the 
asymmetry alon g th e nodal direction with respect to the 
(7r/2,7r/2) point |24l. [ist . and a renormalized bandwidth 
of nearly 1 eV are all faithfully reproduced in the the- 
oretical results. Further, in the lightly doped (x <C 1) 
cuprates, the added holes would immediately populate 
states near kn = ^{tt/2, n/2): this would give rise to nodal 
quasiparticles whose weight grows linearly with x. In a 
rigid-band picture, the quasiparticle dispersion plotted in 
Fig. 3 suggests the formation of small hole pockets in the 
deeply underdoped regime (see the dotted contour), as 
inferred from recent magnetoresistance oscillation mea- 
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surements [29j, the appearance of a holelike cylindrical 
"Fermi surface" at intermediate dopings, and a change 
to an electronlike FS in the overdoped samples (the full 
lines in Fig. 3). Interestingly, a nearly rigid shift of the 
chemical potential was clearly observed with doping in 
several cuprates [s^ 31 1. The renormalized FS which 



would result at intermediate dopings (the dashed curve in 
Fig. 3) shows all the characteristics of that measured ex- 
perimentally (3, [2HI : correlations induce pronounced flat- 
tening of the dispersion near the antinodal points, while 
electronic states in the nodal region remain dispersive. 
That the ZR-z^ mixing is important is shown by the fact 
that it drastically reduces the dispersive character close 
to fcan = (tt/o, 0), bringing the computed ZR dispersion 
as well as the FS in much closer agreement with the ex- 
periment. The band at higher binding energy, associated 
with the charge transfer from the d322_^2 level to the 
in-plane p-d orbitals, was not resolved experimentally. 
More conclusive elaboration of this aspect requires one 
to explicitly treat the dynamical spectral-weight transfer 
between the ZR and higher excited states: whether this 
can lead to an understanding of "waterfalls" in the dis- 
persion at higher binding energy is out of scope of the 
present work. 

Our theoretical findings have other important implica- 
tions. Very recent electronic Raman scattering measure- 
ments [13] clearly show the "nodal-antinodal" dichotomy 
in cuprates. For small x, the nodal lineshapc exhibits 
a low-energy increase, roughly ~ x, while the antinodal 
lineshape is almost independent of x, fully in agreement 
with our results above. Further, the shape of the derived 
FS implies existence of the "nesting" vectors, Q, joining 
the flat-band parts of the renormalized dispersion, see the 
inset in Fig. 3: these are important when a description 
of the inelastic neutron scattering (INS) measurements 
is attempted Uj. ARPES, Raman, and INS results can 
thus be reconciled with our calculated FS. 

To conclude, quantum chemical MC calculations were 
employed to study the lowest electron-removal states in 
an undoped layered cuprate. The method is sufficiently 
general and readily adaptable for studying the correlated 
band structure of other 3d-oxides of intense current inter- 
est. The ZR physics and farther spin interactions around 
a ZR hole are treated on equal footing in our approach, 
which allows an accurate description of the renormaliza- 
tion effects on the dispersion of the ZR-like quasiparti- 
cle band. FM correlations involving spins at nearest Cu 
sites and AF couplings beyond these nearest neighbors 
are both treated consistently. The mixing with "triplet" 
d^z2_^2 hole states at higher binding energy is also in- 
vestigated. We show that all these effects are important 
for describing the dispersion of the lowest (A^— 1) band. 
In contrast, the FM d-d correlations induced by the ZR 
2p hole and the ZR-z^ mixing are missing in standard 
t-J or single-band Hubbard models. The agreement be- 
tween our theoretical results and the available ARPES 



data is remarkable. We are able to reproduce all the im- 
portant details of the experimental spectrum. Extension 
of our ab intio approach to tackle dynamical responses is 
an outstanding open problem, which we leave for future 
consideration. 

We thank S. Nishimoto and P. Fulde for insightful dis- 
cussions. 
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